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Abstract — The acquisition stage in GNSS receivers determines 
Doppler shifts and code phases of visible satellites. Acquisition 
is thus a search in two continuous dimensions, where the digital 
algorithms require a partitioning of the search space into cells. 

We present analytic expressions for the acquisition perfor- 
mance depending on the partitioning of the Doppler frequency 
domain. In particular, the impact of the number and width of 
Doppler bins is analyzed. The presented results are verified by 
simulations. 

Index Terms — GNSS, acquisition, Doppler bin width, receiver 
operating characteristics 



I. Introduction 

In Global Navigation Satellite Systems (GNSS) every satel- 
lite is transmitting a particular pseudo-random noise (PRN) 
code, which is known at the receiver Satellites are acquired 
by correlating the received signal and local code signals and 
comparing the results against a threshold. In practice, the local 
replica of the transmitted code signal differs from the received 
code signal by a code phase shift (i.e., time lag) and a Doppler 
shift. Both have to be determined simultaneously in a two- 
dimensional search. The results of this search, which is usually 
called acquisition, are required for presetting subsequent stages 
of the GNSS receiver 

For this two-dimensional search, the continuous time- 
frequency uncertainty region is divided into cells, each cor- 
responding to a particular Doppler frequency and a partic- 
ular code phase. Typically, the number of considered code 
phases is predetermined by the sampling rate and optional 
decimation/interpolation methods, whereas the width (and, 
thus, the number) of Doppler bins is only limited by the 
effective bandwidth of subsequent signal processing stages [ 1 ] . 
In many civil GNSS receivers exploiting the GPS LI C/A 
code fT\, the integration period of the correlator is set to the 
code period of 1 ms. The corresponding Doppler bin widths 
range from 500 to 667 Hz (see Il3l-||5] as well as (|6l and the 
references therein). Unsurprisingly, the choice of the Doppler 
bin width strongly influences the acquisition performance: Not 
only that the Doppler bin width is inversely proportional to the 
number of cells to be searched, it also strongly influences the 
probability of signal detection. Aside from that, the probability 
of false positive detections at a Doppler bin adjacent to the 
correct one increases for small Doppler bins. 



Non-coherent acquisition methods take the squared magni- 
tude of the correlation coefficients as a decision metric to over- 
come unknown carrier phases and possible data modulation. 
For these methods, the computation of the receiver operating 
characteristics is a well-investigated field of research. The liter- 
ature provides detection and false alarm probabilities for single 
cells ||5|, a serial search over all cells with threshold compar- 
ison Q, El, a maximum search [9|, [10] and combinations 
thereof [11]. In [12], a comparison of the abovementioned 
techniques is provided for an LI GPS receiver Detection 
probabilities for an L5 GPS receiver with different algorithms 
combining data and pilot signals are considered in ifTSJI . from 
which the comprehensive signal model was largely adopted in 
this work. In [12], both the number of Doppler bins and side 
lobes resulting from adjacent Doppler bins are considered; the 
latter are only obtained by means of simulations. Analytical re- 
sults for the effect of residual Doppler shifts on the acquisition 
performance have been presented in llT4l . There, also the effect 
of a single adjacent Doppler bin containing significant energy 
was analyzed, however, a detection in either of these bins was 
considered correct and detection performance was evaluated 
only numerically. All these works, however, are lacking an 
analysis of the influence of the Doppler bin width on detection 
performance in terms of closed-form expressions. 

In this work, we fill this gap by deriving expressions for 
cell detection probabilities as a function of the Doppler bin 
width. These cell probabilities are then used to compute 
global detection and false alarm probabilities, which further 
depend on the number of Doppler bins. With the help of 
this theoretical framework, a proper analysis and, maybe even 
more importantly, a performance-oriented design of GNSS 
acquisition stages is possible. Moreover, while the focus of 
this work is on GNSS receivers, the results can be applied to 
other CDMA systems affected by large Doppler shifts. 

The remainder of this article is organized as follows: In 
Section |II] the signal model is introduced, while Section |III] 
gives a detailed analysis of the acquisition process. The main 
contribution of this work is concentrated in Sections HVlandrVl 
The former is devoted to deriving global detection and false 
alarm probabilities for generalized cell probabilities, while 
in the latter the influence of the Doppler bin width on cell 
detection probabilities is discussed. The analytic results are 
finally verified by extensive simulations in Section IVII 
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II. Signal Model 

After front-end filtering, downconversion to the intermediate 
frequency (IF), and A/D conversion, the signal received from 



single satellite can be represented as lfT2l 



riF[n\ — V2C y[n\ cos 



9iF + 0D)n-d 



V[n\ 
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where C is the signal power, 6jp = '^^y'' and 
are the sampled angular equivalents of the intermediate and 
Doppler frequencies fip and fn, respectively, and i9 is a 
phase shift introduced by the transmission and the noncoherent 
downconversion to IF. If the ideal front-end filter has a band- 
width equal to ^- with fs being the sampling frequency, the 
noise signal ri\n] is assumed to be Gaussian with variance ifTSl 
pp. 556, Prop. 25.15.2] 



a.,-, 



(2) 



and with an autocorrelation function of ifTSI 



^?),?)[™] = E{i][n]r][n — m]} = cr,'J5[m]. (3) 

In these equations, ^ is the two-sided noise power spectral 
density. The carrier is modulated by 

y[n] = d[n]c[n] (4) 

where d[n] is the data message and c[n] is the binary 
PRN code, i.e., c[n] — ±1. For the sake of simplicity, this 
work only considers the GPS LI C/A code with a number 
of chips per code period Nq = 1023 and a code period 
Tper = 1 ms |5|. Furthermore, it is assumed that no data is 
modulated on the PRN code, i.e., d[n] — 1. This assumption 
is unproblematic for the GPS LI C/A codes, since the 20 ms 
duration of a data bit is significantly larger than the code 
period. If acquisition is performed twice for two consecutive 
C/A code periods and by taking the result with the stronger 
correlation peak, it can be guaranteed that no bit transition 
occurs within the considered code period (l3|. For modern 
codes, where the code period is equal to the bit duration 
(e.g., the GPS L2 CM code) the effects of bit transitions can 
be mitigated using non-coherent integration lfT6l or aided 
acquisition ifTTl . 



III. Acquisition System 

As stated in the introduction, the continuous time-frequency 
uncertainty region has to be partitioned into cells to make 
acquisition tractable. Let us, for the remainder of this work, 
assume that the partitioning of the frequency domain is uni- 
form, and that each of the K resulting Doppler bins has a 
width W, as shown in Fig. [T] Given a maximum expected 
Doppler frequency ±fD,ma.x, the number of bins, K, is given 
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In the acquisition process illustrated in Fig. |2] the received 
signal rjF [n] is first downconverted using an expected Doppler 
frequency 6^,. The obtained signal rgln] can be described 
utihzing ([T]i by 



rB[n]^riF[ny'-'"^+'"^>'' (6) 
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Fig. 2. Acquisition of a signal with unknown Doppler frequency and code 
phase. 



with Mf, 



'^. — Qd and YiOj, — 9f^ + 6d- The noise signal 
fjln] is a zero-mean circular-symmetric complex Gaussian 

2 

(ZMCSCG) signal with variances ^ for real and imaginary 
parts. After downconversion the signal is multiplied with the 
spreading code using an expected code phase m. Thus, 



rB[n\c[n - 
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where c[n — rh] is the code c[n] circularly shifted by to. The 
decision metric X[Th,k] is obtained by averaging the signal 
over one code period Tper, 



1 ^-^ 
X[m,k] = ^ X! ^t*^] 

n=0 



(10) 



where the number of samples within one code period is given 
by 

N = Tperfs- (11) 

In this operation, the high-frequency term ^Ji'^^iF+^^k)'^-o^ 
vanishes and the bounded sum over the low-frequency term 
can be represented by a Dirichlet kernel: 
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Following the reasoning in ifTSll . the influences of time lags 
and Doppler frequencies can be separated on average. This 
argument is supported by extensive simulations showing that 
the error resulting from this approximation is well below 
maximum side lobe levels (-21 dB according to |5|) and can 
thus be neglected. We therefore get for the decision metric: 



X[jn, k\ — e 
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where Ry^clm] is the correlation function between y[n] and 
the local code c[n\ evaluated at lag rh. The noise signal n[rh] 
is the average of N independent ZMCSCG samples, thus the 
variances of the real and imaginary parts reduce with (fTTT i to 
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Note that the spreading code c[n] has to be upsampled to 
the sampling rate fs prior to correlation. Depending on the 
implementation of the correlation (matched filter, parallel code 
phase search |3|, etc.) samples of the decision metric adjacent 



-t per 



O 

Q 















m 












Ts 


































































































































































































1 





















































































































































k+1 

k 

k-1 



Vd, 



N code phases m 

Fig. 1. Partitioning of tlie continuous time-frequency uncertainty region. Tlie coirect code pliase is denoted by m, the coirect Doppler bin by k. The plus-sign 
indicates the true Doppler frequency 9n, the dot shows the Doppler estimate d^^ minimizing the residual Doppler difference A6^. Note that adjacent Doppler 
bins can contain significant signal energy (indicated by different shading; cf. Sections llVI and fvl. 



in the code phase domain are not necessarily statistically 
independent: The correlation function can, e.g., be computed 
for all lags m and a particular Doppler estimate k based on 
the same set of A^ input samples. In this case the correlation 
sum runs over differently weighted moving average filtered 
input samples, where the weights are ±1 and the filter order 
is identical to the upsampling factor -^ lfl9l . This dependence 
can be exploited to significantly reduce the computational 
complexity of the correlation process by means of averaging 
correlation ll20l . 11211 . The result of this averaging correlation 
is not only a decimation of the correlation function from 
sampling rate to chip rate (i.e., N — Nc), but also yields 
statistical independence of samples adjacent in the code phase 
domain. 

The decision is finally based on the squared magnitude of 
the decision metric X[Th,k], 



\X\m,k] 
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which follows for given m and dj^ a non-central x^ -distribution 
with two degrees of freedom and the non-centrality parameter 

.hi{X[m,k]}\ E^^{X[m,k]}\ 
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Note that the squared means have to be normalized by the 
corresponding variances, since the x^ -distribution is defined 
as the sum of squares of Gaussian random variables with unit 
variance 1,22, pp. 940]. With 



A% = 27r 






(18) 



and due to the fact that the maximum Doppler difference A/^, 
is in the order of a few kHz, even for N = Nc = 1023 the 
Dirichlet kernel can be well approximated by a sine kernel. 
Using sinc(a;) = '^'"^^^^ this finally leads to 



C 
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IV. Detection Probabilities 

We define a global detection as the event that the cell 
selected by the employed search strategy is the correct cell, 
i.e., the one with the correct code phase rh — ra and with the 
correct Doppler index 



arg min { | i 
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h\] 
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If a cell is chosen in the absence of a signal, we will call this 
event a global false alarm. If the search algorithm chooses no 
cell at all or the wrong cell in the presence of a signal, neither 
a false alarm nor a detection occurs. We will limit ourselves to 
search strategies employing threshold comparison, i.e., a cell 
detection or cell false alarm is triggered whenever the decision 
metric |X[77i, k] p for this cell exceeds a certain threshold /3. A 
decision based on the ratio between the largest and the second 
largest value of |X[m, A;]p of a subset of cells was suggested 
by 1231 , Il24ll after the introduction of this ratio as a reliability 
measure in ifTTl . However, the performance of this method 
has been analyzed just recently ||25l , and it was shown that a 
decision based on threshold comparison outperforms the ratio 
detector |26|. 

On one hand, as it can be seen from ( fT9] l, the non-centrality 
parameter L^.^ f, of the x^ -distribution is maximized for the 
correct code phase m and the correct Doppler bin k. On the 
other hand, whenever the desired satellite PRN code sequence 
is not contained in the received signal, whenever the difference 
between the actual and the estimated Doppler frequencies is 
too large, or whenever the code phase is not correct (m ^ m), 
the non-centrality parameter L^ fe ~ 0, neglecting side lobe 
and cross-correlation levels. Assuming L^ fc = immediately 
translates to the fact that the x'^ -distribution changes from a 
non-central to a central distribution. 

In this Section, we will derive general relations between 
the cell and the global detection and false alarm probabilities. 
In particular, we consider not only the number K of Doppler 
bins, but also take into account that cells with correct code 
phases, but wrong Doppler indices, may have a non-centrality 
parameter L^j, > 0. 



A. Cell Detection Probabilities 

Let us define the cell false alarm as the decision metric 
|X[m, fc]p exceeding a certain threshold /3 given that L^~_^ f, = 
0. Then, the cell false alarm probability becomes 

F/a(/3) = Prob {\X[m, k]\^ > (3 \ L^j^ = o) (21) 

which for a central ^^-distribution is equivalent to ifTOl 



PfaiP) 



(22) 



Conversely, whenever the decision metric exceeds the thresh- 
old for a non-central ^^-distribution {L^ f, > 0), we will call 
this event a cell detection. Thus, the cell detection probability 
is 



PdetiP.L. 



m,fc]p>/3|i. ,>0 




(23) 



where Qi (•, •) is the Marcum Q-function ||27l, ||28l. Note that 
in fact Pfa{f3) is a special case of Pdet{(3, L^ ^), namely for 



L 



m , k 



0. 



B. Global Detection Probabilities - Naive Assumption 

Let us now assume that the acquisition is implemented as a 
serial search over the two-dimensional uncertainty region com- 
prised of NK cells, and the search is stopped when |X[to, fc]p 
exceeds the threshold /3 for the first time. Furthermore, it 
is assumed that there is only one cell containing significant 
signal energy, i.e., there is only one (to, fc)-pair for which the 
decision metric |X[to, fc]p is non-centrally x^ -distributed - 
namely the pair (to, k). We will call the event of a threshold 
crossing a global false alarm whenever the desired PRN code 
is not contained in the received signal. A global detection, as 
already mentioned, denotes the event when the first threshold 
crossing in the serial search occurs at the correct code phase 
TO, and the correct Doppler bin k. Following [IZ], the global 
false alarm probability Pfa{P) therefore calculates to 



PpaW) = I - a - Pfai(3)Y 



(24) 



whereas the global detection probability Pdet{I3) can be 
calculated as 

sNK 



P. 



DET 



m = 



1 i-ji-Pf^m 
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-Pdetil3,L^,k)- (25) 



Pfail3) 

For very small values of PfaiP) the above equations can be 
approximated by PfaW « NKPfa{(i) and PoETifi) ~ 

Pdet{P) ma. 

The assumption of a single signal cell (i.e., a cell for which 
^rh k ^ ^) ^^ clearly a strong one, since the side lobe and 
cross-correlation levels of the correlation function, a non-zero 
width of the correlation main lobe for N > Nc, and signal 
energy in Doppler bins adjacent to the correct one will lead 
to L^ fe > for more than one cell. Effects of side lobes and 
cross-correlations in the code phase search can be neglected 
in medium SNR levels or mitigated by appropriate threshold 
settings, and the effects due to the correlation main lobe can be 



reduced by means of averaging correlation lfT9l - ll2n . Signal 
energy in adjacent Doppler bins, however, not only affects 
global detection probabilities, but also depends on the width 
of the Doppler bins. A proper analysis of this influence can 
not be found in the literature. 

C. Global Detection Probabilities - Refined Model 

In accordance with what has been said, let L„i fe be the 
non-centrality parameter of the correct cell, and let £,,,.^±1, 
Lm.k±2, etc. denote the non-centrality parameters of cells with 
the correct code phase to, but for the first, second, etc. adjacent 
Doppler bins. For all other code phases let L^ fc = (in the 
absence of the desired PRN code let L^ ^ = for all cells). In 
other words, each Doppler bin k contains at most one signal 
cell at the correct code phase to (see Fig. ^ with an expected 
non-centrality parameter L^ f,. Again, acquisition is assumed 
to take place as a serial search with threshold comparison. It 
is important to note that increased detection probabilities for 
wrong Doppler bins are adverse in terms of global detection 
probabilities. 

Since there is now more than a single signal cell, the 
direction of the serial search has an influence on the global 
detection probability. In other words, searching all code phases 
for each Doppler bin and searching all Doppler bins for each 
code phase leads to different performance results. Note that 
both methods can be efficiently implemented using the FFT: 
In this case, the former option is called parallel code phase 
search, while the latter is often referred to as parallel frequency 
search ID. Following the reasoning in lfT2l . the probability of 
detection for a search over all code phases for each Doppler 
bin can be calculated as 
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Where P/,(/3) = 1 - P/,(/3) and Pdet(/3, i„,fc) = 1 - 
Pciet{l3,L^ ^). The probability of detection for a search over 
all Doppler bins for each code phase calculates to 
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Although theoretically different, for small cell false alarm 
probabilities Pfa{P) we can write Py^ sa 1 — MPfaiP) w 1 
and thus obtain as an approximation for both search directions: 
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(28) 



Depending on the acquisition strategy, this approximation can 
be shown to hold for a wide range of SNR. This is due to 



setting the threshold to obtain a constant false alarm rate ||3], 
which can be done efficiently by continuously measuring the 
noise floor |26| . 

It is worth mentioning that the global false alarm probability 
is the same for the naive and the refined model, i.e., it is 
always given by ( l24b . Further, when the naive assumption 
holds, i.e., L^ ^ = foi" rh ^ m and k ^ k, both ( |26] l 
and (l27b reduce to dZST l. while the same does not hold for 
the approximation ( |28] |. Finally, both the naive and the refined 
model assume that the correct cell is uniformly distributed 
over the two-dimensional search space. While this assumption 
can be justified for the code phase domain, depending on 
the approximate time and position of the receiver, an esti- 
mated Doppler frequency can be computed, from which the 
search should be initiated [l]. Since the Doppler frequency 
is likely close to its estimate if the estimation process was 
successful, the uniformity assumption is too restrictive and 
the obtained results underestimate the detection performance 
of the receiver. 

V. Influence of Doppler Bin Width on the Decision 
Metric 

We have defined the non-centrality parameter of the x^- 
distribution modeling the decision metric |X[to, fc]p in ( fT9] l, 
which shows that it is proportional to the squared sine of the 
Doppler difference A/j. = /^ — /d, as well as to the squared 
correlation function R^ ^[rh]. Using proper decimation meth- 
ods lfT9]| - ll2Tl and neglecting side lobes and cross-correlation 
levels we assume that there exists only one single code phase 
m for which the correlation function is non-zero, and that for 
this phase we have Ry^^lm] = 1. In the literature (e.g., ID, 
HJI), however, even the squared sine is often approximated 
by unity for the correct Doppler bin, and by zero for all other 
bins such that 



^m.k 
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(29) 

- =0 Vfc 7^ fc. (30) 

As Fig.|3]shows, this simplification is too optimistic: For small 
Doppler bins, such as W^ = 200 Hz, the Doppler bin adjacent 
to the correct one (the area between the two leftmost dot- 
markers) has a significant non-centrality parameter Lm.k±i- 
Moreover, for large Doppler bins (e.g., W — 700 Hz) the non- 
centrality parameter L„j ^ of the correct Doppler bin might be 
as low as indicated by the leftmost cross-marker This Section 
is thus devoted to a more in-depth analysis of the influence 
of Doppler bin widths, which will lead to a more realistic 
characterization of the acquisition performance. 

As stated in Section |III] (see Fig. [TJ, in the acquisition 
process the whole Doppler domain is divided into bins of equal 
width W, which are searched in a serial fashion. The center 
frequency of a bin represents the Doppler estimate /j. = df.-^ 
with which the received signal is demodulated. In the Doppler 
bin corresponding to the correct Doppler frequency, the re- 
maining difference A/^ is uniformly distributed within that 
bin, i.e., 

Mk-Ui- — , — ]. (31) 
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W = 200 Hz 
O W = 500 Hz 
X W = 700 Hz 




1 000 1 500 

Doppler shift in Hz 

Fig. 3. Squared sine function for Tper = 1 rns. Markers indicate Doppler 
bin boundaries for designated bin widths W . 



Further, due to the symmetry of the sine, for the k ± ^-th bin 
the residual Doppler difference is distributed according to 



A/fe±i ^U[{21 



w w 



(32) 



with I E N. Since now A/j, is a random variable (RV), also L 
becomes an RV with probabihty density function (PDF) fL{l)- 
The detection probability thus calculates to 



Pdet{l3,L) 



fY\L[v\l)fL{l)dldy (33) 



where /y|l(2/|0 is the PDF of |X[to, fc]p for the cell con- 
taining the desired signal conditioned on the non-centrality 
parameter L. This expression is difficult to compute since 
the PDF of L is not readily available and also the resulting 
integrals might not have convenient closed-form solutions. 
Instead, assuming that the Doppler bins are sufficiently small, 
it is possible to approximate the PDF of Y to be linearily 
dependent on L and we get 



Pdetili.L) 



/y|L(y,E{L})dy, 



(34) 



P 



as shown in the Appendix. Exploiting the method for com- 
puting the expected value of a function of a random variable 
from |29^, pp. 142] and using the uniform distribution (|32] |. it 
is straightforward to compute L 
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Substituting x ~ Aff^Tpe,- we get dx = Tperd{Af^.) and 
2 c rfi'^p' 
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where /; = (2/ - l)-f and /; == (2/ + l)-f denote the lower 
and upper frequency bounds of the k ± ^th Doppler bin. We 



obtain 
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where Si(-) is the sine integral ll22l pp. 231]. For f, = fi = ^ 
(i.e., for the bin k containing the correct Doppler frequency) 
this yields 
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TrVFTper 



(38) 



The sine integral increases monotonously in the interval 
[— TT, tt] and oscillates around the constants -^ for positive 
and around — ^ for negative arguments in a decaying manner. 
Thus, the difference of the sine integrals in ( [37] ) will contribute 
significantly whenever at least one of the arguments falls 
within [— 7r,7r], i.e., for low values of / = |fc — k\. Thus, one 
can expect that Doppler bins close to the correct bin have, on 
average, large non-centrality parameters, while large offsets 
lead to small non-centrality parameters. These mathematical 
considerations are in line with intuition, which suggests that 
the correct Doppler bin contains most of the energy. Moreover, 
since for small values of W more bin boundaries may fall 
in the interval [— 7r,7r], the correct Doppler bin k on average 
contains more energy if the bin is small. For large bins the 
correct Doppler frequency may be far from the bin center, 
compared to the former case. 

While all this reasoning suggests that smaller Doppler bins 
are preferable, another aspect has to be taken into account: 
Given a fixed Doppler search range of ±/D,max, according 
io ^ K bins have to be searched. Since K is inversely 
proportional to the bin width W, smaller widths lead to a 
higher number of bins, which increases the probability of 
false alarms according to (l24l l and, consequently, decreases the 
probability of detection (see d25T l and (|28]l). Moreover, in cases 
where W is small, bins adjacent to the correct Doppler bin may 
contain significant signal energy (large Lm,k±i, Lm.k±2, etc.), 
which can trigger a false detection and thus degrade receiver 
performance. As a consequence, depending on the threshold 
and the SNR, there will be an optimal Doppler bin width W 
maximizing the global detection probability Pdet{P)- 

VI. Simulations and Results 

To verify the analytic results, a series of simulations was 
conducted. To this end, a set of satellite signals was generated. 
For simplicity, it was assumed that just a single satellite (GPS 
LI C/A PRN code 1, Nc = 1023, Tper = 1 ms) was 
visible with random Doppler frequency fo and code phase m 
(see below). The carrier-over-noise spectral density ratio -^ 
was set to 40 dBHz unless stated otherwise. After sampling 
the signal with a high sampling frequency, it was assumed 
that prior to detection the signal was decimated to the code 
chipping rate by means of averaging correlation (cf. Il20l . 
Ell). Thus, N = Nc ^ 1023. This simpHfication does not 
affect the validity of the analysis, since with this method the 
statistical properties of the cells do not change ||T9| . It was 



assumed, however, that during decimation the correct code 
phase m corresponding to Ry^dm] = 1 is preserved. 

The Doppler frequency was assumed to be uniformly 
distributed over the whole Doppler range, i.e., fn ^ 
l^ (-/D,max, /D,max), whcrc /zj^max = 5000 Hz |1|. To make 
both simulation and analytic comparison tractable, only two 
Doppler bins adjacent to the correct bin contained signal 
energy, i.e., Ljn,k±i — for Z = 3, 4, . . . , which leads to 



Pdet{/3,L,„,k±l) - Pfa (/?) VZ = 3,4, . . 



(39) 



This simplification holds well for bin widths greater than 
300 Hz, as shown in Fig. [3] The Doppler bin widths were 
vai-ied within W S {200, 500, 700, 1000}. The signals were 
correlated with PRN codes 1 and 5 for the detection and 
false alarm probabilities, respectively. For each Doppler bin 
width W and each value of -;^ a set of 10^ correlations was 
performed. 

A serial search is implemented: Starting from the first 
Doppler bin, all possible code phases are searched sequentially 
until either the threshold is crossed or until the whole Doppler 
bin is searched. Then, the next Doppler bin is taken into con- 
sideration. If the first threshold crossing occurs at the correct 
code phase in the correct Doppler bin, the signal is assumed 
to be detected, while any threshold crossing in the absence of 
a signal triggers a false alarm. For the analytic results for the 
global detection probability, the accurate expression (|26] | was 
used. 



A. Influence of Doppler Bin Width on Cell Probabilities 
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Fig. 4. Cell Detection Probability Pdet{P, Lm,k) for different Doppler bin 
widths. Simulated (bold markers) and analytic (lines) results are shown for 
the connect Doppler bin. The thick solid line indicates Pdetil^, imax). 

Figs. mis] and |6] show a comparison between the simulated 
and analytic cell detection probabilities for the correct, the 
first, and the second adjacent Doppler bins for different bin 
widths. It can be seen that there is a good match between 
the analytic and the simulated results, except for the case of 
the Doppler bin directly adjacent to the correct one. In this 
particular case, a separate analysis showed that for k ~ k±l 
and for larger the Doppler bin widths, (|34] | is not a good 
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Fig. 5. Cell Detection Probability PdetC/^i im,fc±l) f™ different Doppler 
bin widths. Simulated (bold markers) and analytic (lines) results are shown 
for the Doppler bin adjacent to the correct one. The thick solid line indicates 

Pdet{l3,0) = PfaW). 



Fig. 7. Global False Alarm Probability Ppa{I3) for different Doppler bin 
widths. Simulated (bold markers) and analytic (lines) results are shown. 
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Fig. 6. Cell Detection Probability PdetC/^i im,fc±2) f™ different Doppler 
bin widths. Simulated (bold markers) and analytic (lines) results are shown 
for the Doppler bin second adjacent to the correct one. The thick solid line 



indicates Pdet(/3,0) = P/,(/3). 



approximation of (|33| |. which consequently leads to large devi- 
ations. This is related to the fact that the linear approximation 
of the conditional PDF fyiLillil) is not sufficiently accurate 
for these choices of parameters (see Appendix). Nevertheless, 
it can be seen that smaller Doppler bins lead to increased 
cell detection probabilities, both for the correct and adjacent 
Doppler bins. In addition to that, it is shown (thick Unes) that 
by not considering Doppler bin widths at all, the results would 
be too optimistic, leading to an overestimation of the global 
detection probability (see Section IVI-Bb . 

The cell false alarm probability does not depend on the non- 
centrality parameter L, thus it is not affected by the Doppler 
bin width. 

B. Influence of Doppler Bin Width on Global Probabilities 

Combining these probabilities to global detection and false 
alarm probabilities shows another picture: Here, the effect of 
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Fig. 8. Global Detection Probability Pdet(P) for different Doppler bin 
widths. Simulated (bold markers) and analytic (lines) results are shown. 



a greater number of noise cells becomes apparent, showing 
that smaller Doppler bins do not necessarily lead to improved 
performance. For example. Fig. |7] shows that the global false 
alarm probability increases for smaller Doppler bins, i.e., for 
an increasing number K of bins - this is intuitively understood 
by looking at ( |24] |. Conversely, Fig. [8] shows the probability of 
detecting the correct cell in a serial search which, especially 
for small Doppler bins and low thresholds, suffers from 
high false alarm rates. The additional bend in the curve for 
W = 200 Hz near the maximum is due to significant energy 
of the adjacent Doppler bins, which increases the probability 
of triggering a false alarm. As shown, low thresholds benefit 
from larger Doppler bins (small K, little energy in adjacent 
bins), whereas the opposite is true for larger values of /3. 
There, small Doppler bins lead to a high imfc, i.e., a high 
cell detection probability for the correct cell, whereas false 
alarms are unlikely due to the large threshold. 

As Figs. I2] and |8] show, the analytic results are widely 
validated by the simulations, despite the fact that the detec- 
tion probability for the Doppler bin with index fc ± 1 was 



underestimated by the theoretical approximation. 
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Fig. 9. Global Detection Probability Pdet{I^) f™ different Doppler bin 

widths as a function of -§-■ Analytic results are shown for a fixed cell false 

No •" 

alarm probability of PfaW ~ 10"^ (/3 = 0.0225). Thick lines coiTespond 
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: for all kj^k. 



Fig. |9] shows the global detection probability for a fixed 
cell false alarm probability and different values of -2-. It 
can be seen that for low values of -^ all Doppler bin 
widths perform almost identically with smaller bins slightly 
in favor. Differences only become apparent for medium to 
high values. Fixing the threshold relative to the noise floor has 
the effect of increasing the probability of a false alarm at an 
adjacent Doppler bin for large SNR. The detection probability 
asymptotically approaches the probability that the serial search 
is started at the correct bin, i.e., Pdet{P) — j^- Thus, for high 
SNR and a fixed cell false alarm probability, large Doppler 
bins perform better than smaller ones. 

As one can expect, by neglecting the influence of the 
Doppler bin width on the cell probabilities, i.e., assuming 
that there is just one signal cell in the two-dimensional 
uncertainty region, this behavior cannot be observed. As the 
thick solid lines in Fig.|9]show, the global detection probability 
just depends on the SNR and on the total number of cells, 
NK. Thus, for lai'ge -^, neglecting adjacent Doppler bins 
leads to optimistic results for the global detection probabiUty 
PdetW). 

Fig. [To] depicts the global detection probability as a function 
of the Doppler bin width W, again for a fixed cell false alarm 
probability. It can be seen that small Doppler bins are in favor 
for low values of -^ (see Fig. |9]). If the bins are large, too 
much energy is lost for the correct bin, while the energy from 
adjacent bins is buried under the noise floor anyway. For high 
SNR regions and small W, however, false alarms are likely 
at adjacent Doppler bins. Again, by neglecting the influence 
of W on the cell detection probabilities, too optimistic results 
are obtained. As the thick lines in Fig. [TO] show, one may 
be falsely led to the conclusion that the effect of Doppler 
bin widths is negligible, and that regardless of the SNR the 
detection performance is better for large Doppler bins (i.e., for 
small K). 




- 37 dBHz 
40 dBHz 
43 dBHz 



200 300 400 500 600 700 800 900 1000 

Doppler Bin Width W in Hz 

Fig. 10. Global Detection Probability Pdet{I3) for different -^ as a 
function of Doppler bin widths. Analytic results are shown for a fixed cell false 
alarm probability of P/a(/3) ~ W~^ (/3 = 0.0225). Thick lines correspond 
to Lm,k = imax, L r = for all k ^ k. 



Here, we want to stress that a general statement about the 
optimal width W and the corresponding number K of the 
Doppler bins is not possible. The naive assumption of a single 
signal cell only considers the number K of bins, and thus 
suggests that few large bins outperform many small ones. 
This effect of the Doppler bin count K is outweighed by 
considering the non-centrality parameters of adjacent bins - 
if the bins are too small, false alarms at adjacent Doppler bins 
prohibit high detection probabilities. As it turns out, there is 
a Doppler bin width W optimizing the trade-off. A general 
statement about this optimum is not possible, since it depends 
not only on the SNR (cf. Fig. [TOl ), but also on the chosen cell 
false alarm probabiUty which influences the threshold setting. 



C. Receiver Operating Characteristics 
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Fig. 11. Receiver Operating Characteristics for 

(red bold markers) and analytic (marked lines) results are compared to the 
naive assumption of L^j^ = Lmax and L r. = for all fc 7^ fc (thick 
marked line) 



Fig. [TT] finally plots the global detection probability as a 
function of the global false alarm probability. Again, it can 
be seen that analysis and simulation match quite well. As it 
was expected, neither very small nor very large Doppler bins 
perform well: In the former case, the false alarm probability 
is increased, while in the latter the detection probability is 
decreased. On the other hand, medium Doppler bin widths 
in the order of 500 to 700 Hz turn out to perform optimally 
in terms of receiver operating characteristics, at least for this 
particular value of -^. 

By neglecting the influence of Doppler bin widths, one again 
can see that the obtained results are overly optimistic. Looking 
at the thick lines in Fig. (TT] which consider the number of 
bins, K, but not the influence on the non-centrality parameters 
L^ y one is tempted to conclude that large Doppler bins out- 
perform smaller ones. This result, however, is based entirely 
on the assumption that there is just one signal cell, and that 
the energy contained in this cell is independent of the Doppler 
bin width. 

As Fig. [TT] shows, despite model inaccuracies pointed out 
in Section [VI-AI and in the Appendix, our theoretical frame- 
work matches the simulations quite well. In addition to that, 
compared to the naive assumption of a single signal cell, our 
model leads to dramatically improved estimates of receiver 
operating characteristics. 

VII. Conclusion 

In this article, the influence of Doppler bin widths on GNSS 
acquisition performance is analyzed. Analytic expressions, 
linking the Doppler bin width to the detection probabilities, 
are derived and evaluated. These expressions replace the 
conventionally used detection and false alarm probabilities, 
which assume that the search region is populated by noise- 
only cells except for a single signal cell. This assumption, 
specifically, is shown to be overly optimistic. 

Three different effects of the Doppler bin width are consid- 
ered: First, the number of bins influencing the total number 
of cells over which a serial search is conducted. Here it is 
shown that a large number of bins increases the false alarm 
probability and decreases the detection probability. Second, 
the influence of the bin width on the detection probability of 
the correct cell was analyzed. It turned out that smaller widths 
improve detection performance significantly. And finally, the 
influence of the bin widths on the false alarm probability, 
especially for cells in the close vicinity of the correct cell, 
where significant signal energy can leak to adjacent Doppler 
bins was discussed. It was shown that smaller bin widths 
increase the probability of stopping the search at a wrong 
Doppler bin, and thus have adverse effects on the detection 
probability. As a consequence, extreme cases of very small or 
very large Doppler bins should be avoided when designing the 
acquisition stage. All these analyses are validated by means 
of simulations. 



that 



fY\L{y,l) - [kl + d)fY[y.^{L)). 

Substituting this into ( l33b leads to 



PdetiP.L) 



P J-oo 

oo />oo 



fY\LiyJ)fLil)dldy 



(40) 



(41) 



{kl + d)fY{y,E{L})fL{l)dldi42) 



fY{y,E{L}) / (kl + d)fL{l)dldiA3) 



(44) 



(45) 



/y(y,E{L})(A:E{L} + d)d2; 



fY\L{y,E{L})dy, 



which is ( l34b . However, in this equation one needs to make 
sure that /y(T/,E{_L}) depends on the expected value of L 
such that 



fY{y,E{L}){kE{L} + d)dy = l 



(46) 



for all E {L}. In Fig. [12] this linearity assumption is illus- 
trated. The range between the rightmost marker and the value 
Lm,k±i/Lniax represents the correct Doppler bin. The range 
between the two rightmost (leftmost) markers indicate the first 
(second) adjacent Doppler bin (fc ± 1 and /c ± 2, respectively; 
compare to Fig. [3]l. As it can be seen, both the range between 
the two leftmost markers and the range between the rightmost 
marker and the boundary can be well approximated by a 
line. Thus, ( |34| | approximates (i33] l well and the corresponding 
detection probabilities are matching simulations (see Figs. |4] 
and O. The range between the two leftmost markers, on the 
other hand, is only approximately linear for W — 200 Hz. 
For bin widths of W = 500 Hz and W = 700 Hz the 
linearity assumption clearly does not hold. As a consequence, 
the detection probability obtained for the adjacent Doppler bin 
is vaUd only for W = 200 Hz (see Fig. O. 
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Appendix 

In this appendix, we show that under a linearity assump- 
tion ( l33T l can be approximated by ( l34l i. Thus, let us first assume 



Fig. 12. Conditional probability /yiLCj/iO evaluated for Tper = 1 nis, 
■^ = 40 dBHz, and y = 0.0225. Markers indicate Doppler bin boundaries 
for designated bin widths W . 
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